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Abstract 

Newton's problem of minimal resistance is one of the first problems of 
optimal control: it was proposed, and its solution given, by Isaac Newton 
in his masterful Pnnctpia Mathematica, in 1686. The problem consists 
of determining, in dimension three, the shape of an axis-symmetric body, 
with assigned radius and height, which offers minimum resistance when 
it is moving in a resistant medium. The problem has a very rich history 
and is well documented in the literature. Of course, at first glance, one 
suspects that the two dimensional case should be well known. Neverthe- 
less, we have looked into numerous references and ask at least as many 
experts on the problem, and we have not been able to identify a single 
source. Solution was always plausible to everyone who thought about the 
problem, and writing it down was always thought not to be worthwhile. 
Here we show that this is not the case: the two-dimensional problem is 
more rich than the classical one, being, in some sense, more interesting. 
Novelties include: (i) while in the classical three-dimensional problem only 
the restricted case makes sense (without restriction on the monotonicity 
of admissible functions the problem doesn't admit a local minimum), we 
prove that in dimension two the unrestricted problem is also well-posed 
when the ratio height versus radius of base is greater than a given quan- 
tity; (ii) while in three dimensions the (restricted) problem has a unique 
solution, we show that in the restricted two-dimensional problem the min- 
imizer is not always unique - when the height of the body is less or equal 
than its base radius, there exists infinitely many minimizing functions. 
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1 Introduction 



Newton's aerodynamical problem, in dimension three, is a classic problem (see 
e.g. m3IU). It consists in joining two given points of the plane by a curve's arc 
that, while turning around a given axis, generate the body of revolution offering 
the least resistance when moving in a fluid in the direction of the axis. Newton 
has considered several hypotheses: that the body moves with constant veloc- 
ity, and without rotation, on a very rare and homogeneous medium of particles 
which are all equal; that the axis-symmetric body is inscribed in a cylinder of 
height H and radius r; that the particles of the medium are infinitesimally small 
and immovable (there exists no temperature motion of particles) ; that collisions 
of the particles with the body are absolutely elastic. Newton has indicated in the 
Mathematical principles of natural philosophy the correct solution to his prob- 
lem. He has not explained, however: how such solution can be obtained; how 
the problem is formulated in the language of mathematics. This has been the 
work of many mathematicians since Newton's time (see e.g. Exten- 
sions of Newton's problem is a topic of current intensive research, with many 
questions remaining open challenging problems. Recent results, obtained by 
relaxing Newton's hypotheses, include: non-symmetric bodies py; one-collision 
non-convex bodies |S]; collisions with friction |S]; multiple collisions allowed |10|: 
temperature noise of particles 1X3 • Here we are interested in the classical 
problem, under the classical hypotheses considered by Newton. Our main ob- 
jective is to study the apparently more simpler Newton's problem of minimal 
resistance for a two-dimensional body moving with constant velocity in a ho- 
mogeneous rarefied medium of particles. The first work on a two-dimensional 
Newton- type problem seems to be where the authors study the problem in 
a chaotically moving media of particles (in the classical problem particles are 
immovable). The results in fTTj were later generalized to dimension three |12| . 
This paper is motivated by the results in |12| : when one considers temperature 
motion of particles, the three-dimensional problem admits only two types of 
solutions; while the two-dimensional case is more rich, showing solutions of five 
distinct types. Here we prove that in the classical framework, with an immovable 
media of particles, also the two-dimensional case is more rich: in certain cases 
of input of data (height H and radius r of the body) the problem is well-posed 
(admit a local minima) without imposing the restriction y(x) > on the admis- 
sible curves y(-). This is different from the three-dimensional classical problem 
or the problem in higher-dimensions, where the restriction y{x) > is always 
necessary for the problem to make sense: without it there exists no strong and 
no weak local minimum for Newton's problem of minimal resistance (see e.g. 
[7|E]). We show that for H > ^r the function y{x) — —x is a local minimum 
for the unrestricted Newton's problem of minimal resistance in dimension two. 
In the restricted case, while in dimension three (or higher-orders) the problem 
has always a unique solution, we prove that infinitely many different minimizers 
appear in dimension two for r > H. These simple facts seem to be new in the 
literature, and never noticed before. 
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2 Restricted and unrestricted problems 



In the classical three dimensional Newton's problem of minimal aerodynamical 
resistance, the resistance force is given by R [{/(■)] = J Q 1+ ^ x yi dx. Minimization 
of this functional is a typical problem of the calculus of variations. Most part 
of the old literature wrongly assume the classical Newton's problem to be "one 
of the first applications of the calculus of variations" . The truth, as Legendre 
first noticed in 1788 (see [2]), is that some restrictions on the derivatives of 
admissible trajectories must be imposed: y(x) >0,xE [0, r]. The restriction is 
crucial, because without it there exists no solution, and the problem suffers from 
Perron's paradox [181 §10]: since the a priori assumption that a solution exists 
is not fulfilled, does not make any sense to try to find it by applying necessary 
optimality conditions. It turns out that, with the necessary restriction, the 
problem is better considered as an optimal control one (see ^3 p. 67] and (17)1. 
Correct formulation of Newton's problem of minimal resistance in dimension 
three is (cf. e.g. f7\\Hjl): 

r x 

J l + u{xY 
y(x) = u(x) , u(x) > , 
2/(0) = 0, y(r) = H, H>0, 

where we minimize the resistance 1Z in the class of continuous functions y : 
[0, r] — > K with piecewise continuous derivative. Here we consider Newton's 
problem of minimal resistance in dimension two (see |17p: 

R[u{ - )] = 1 IT^" min ' 

y(x) = u(x) , u(x) G Q , ^ 
y(0) = 0, y(r) = H, H>0. 

We consider two cases: (i) unrestricted problem, where no restriction on the 
admissible trajectories y(-) other than the boundary conditions y(0) = 0, y(r) = 
H is considered (f2 = R); (ii) restricted problem, where the admissible functions 
must satisfy the restriction y(x) >0,i£ [0, r] (£1 = Mg"). While for the classical 
three-dimensional problem only the restricted problem admits a minimizer, we 
prove in S0|that the two-dimensional problem Q is more rich: the unrestricted 
case also admits a local minimizer when the given height H of the body is 
big enough. In we study the restricted problem. Also in the restricted 
case the two-dimensional problem is more interesting: if r > H, then infinitely 
many different minimizers are possible, while in the classical three-dimensional 
problem the minimizer is always unique. 
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3 General results for both problems 

The central result of optimal control theory is the Pontryagin Maximum Prin- 
ciple |13j , which gives a generalization of the classical necessary optimality con- 
ditions of the calculus of variations. The following results are valid for both 
restricted and unrestricted problems: respectively SI = Rq and ft = R in QJ. 

Theorem 1 (Pontryagin Maximum Principle for Q). If (?/(•), it(-)) is a 
minimizer of problem QJ, then there exists a non-zero pair (tpo,tp(-)), where 
ipo < is a constant and ip(-) £ PC 1 ([0, r]; M), such that the following condi- 
tions are satisfied for almost all x in [0, r]: 

(i) the Hamiltonian system 

= ^(u(x), ipo, ip(x)) (control equation y = u) , 
= — jj^(u(x), ipo,ip(x)) (adjoint system ip = 0) ; 

(ii) the maximality condition 

H(u(x), tp ,ip(x)) = maxH(u, ip , tp(x)) ; (2) 

where the Hamiltonian Ji is defined by 

H{u,ip ,ip) = ip - 1 2 +ipu. (3) 
l + u 

The adjoint system asserts that ip(x) = c, with c a constant. From the 
maximality condition it follows that ipo ^ (there are no abnormal extremals 
for problem QJ). 

Proposition 2. All the Pontryagin extremals (y(-),u(-),ipo,ip(-)) of problem 
IjTjl are normal extremals (ipo ^ 0), with ip(-) a negative constant: ip(x) = —A, 
A > 0, x E [0,r]. 

Proof. The Hamiltonian H for problem QJ, H(u,ipo,ip) — ipo x ■ M a + ^u, does 
not depend on y. Therefore, by the adjoint system we conclude that 

ip(x) = (tt(i), Vo, ^(a)) = > 

that is, ip(x) = c, c a constant, for all x G [0, r]. If c = 0, then ipo < 
(because one can not have both ipo and ip zero) and the maximality condition 
J2J simplifies to 

= max J ^ } U) 

1 + u 2 (x) uen \ 1 + u 2 J 

From (0J we conclude that the maximum is not achieved (it — » oo). Therefore 
c ^ 0. Similarly, for c > the maximum 

Tr — + cu{x) = max <^ = + c 

1 + u 2 (x) «en ( 1 + u 2 
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does not exist, and we conclude that c < 0. It remains to prove that tpq ^ 0. 
Let us assume ipo = 0. Then the maximality condition reads 

cu(x) = max{cu} , c < . (5) 

For Q = R the maximum does not exist, and we conclude ipo ^ 0. For = 
(JSJl imply u(x) = and y{x) = w, w a constant (y(x) — u(x)). This is not 
possible, given the boundary conditions y(0) = and y(r) = H with H > 0. 
Therefore V^o ; there exists no abnormal Pontryagin extremals. □ 

Remark 3. If (y(-),u(-),ipQ,'ip(-)) is an extremal, then {y(-),u(-),jipo,'yip(-)) is 
also a Pontryagin extremal, for all 7 > 0. Therefore one can fix, without loss 
of generality, ipQ = — 1 . 

From Proposition |3 and Remark |3| it follows that the Hamiltonian J2J takes 
the form 

H (u) = l — - Xu , A > . (6) 

1 + u z 

It is not easy to prove the existence of a solution for problem Q with classical 
arguments. We will use a different approach. We will show, following [TJ\, that 
for problem Q the Pontryagin extremals are absolute minimizers. This means 
that to solve problem Q it is enough to identify its Pontryagin extremals. 

Theorem 4. Pontryagin extremals for problem Q are absolute minimizers. 

Proof. Let u(-) be a Pontryagin extremal control for problem (Q. We want to 
prove that 

1 , T 1 

-ax > / „ , ; ax 



/ l + u 2 (x) ~ J l + u 2 (x) 

for any admissible control u(-). Given ijfJJl, we conclude from the maximality 
condition J5J that 

' Xu(x) > — \— - \u(x) (7) 



1 + u 2 (x) 1 + u 2 (x) 

for all u(-) G PC ([0, r], O). Having in mind that all the admissible processes 
(y(-),u(-)) of satisfy 



u{x)dx = \ y(x)dx = y(r) — t/(0) = H , 
Jo 

we only need to integrate Q to conclude that u(-) is an absolute control mini- 
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mizer: 



ii 



1 - 


f u 2 {x) 




1 


H 


- u 2 (x) 




1 


H 


- u 2 (x) 




1 


H 


- u 2 (x) 



1 

o V l + u 2 (x) 



Xu(x) ) dx > j ( — - — j — J7~\ ~ Xu(x) ) dx 



I - — ~-2i„\ dx + A / u{x)dx < j 2~T~\d x ^ / u { x )dx 

= / ^r^cte + XH < I - ^--dir + XH 



o l + u 2 (x) 



□ 

Roughly speaking, Theorem 0] reduces the infinite dimension optimization 
problem Q to the study of a one-dimension maximization problem: 

max H (u) = max j — - — Xu 1 , A > . (8) 

«es! «en [ 1 + u 2 J 

4 Unrestricted problem 

The following standard result of calculus (see e.g. will be used in the sequel. 

Theorem 5. Let n > 2 and £1 C W. be an open set. If f : — > K is n — 1 
times differentiable on Q and n times differentiate at some point ae!l where 
f (k) {a) =0fork = 0,...,n-l and / (n) (a) / 0, then: 

• either n is even, and /(•) /ias an extremum at a, that is a maximum in 
case f( n '(a) < and a minimum in case /' n ^(a) > 0; 

• or n is odd, and /(•) does not attain a local extremum at a. 

We are considering now the unrestricted two-dimensional Newton's problem 
of minimal resistance, that is, = R in 0). A necessary (sufficient) condition for 
u to be a local maximizer for problem JSJ is given by Ti' (u) — and H" (u) < 
(H" (u) < 0), where 

2?; 

W (u) = =^>-A, 

(l + « 2 ) 2 

H"(.) = -2. 3m2 - 1 



(1+M 2 ) 3 ' 

From the first order condition (maximality condition J2Jl) it follows that 

u(s) A y(x) A 

(l + u 2 (a;)) 2 2 {l + y 2 (x)f 2 

Using the boundary conditions y(Q) — and y(r) = i/, we conclude that y(x) — 
—x(u — —)isa. local candidate for the solution of the unrestricted problem 
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= fri+gpy* )• However, by Theorem^ we conclude that such u is a maximizer 
only when H > -y^f*. For < the value u = — corresponds to a local 
minimizer of H (u) since K" > 0; for H — function H (u) has neither local 
maximum nor minimum since TL" (j^f^ =0 and H'" (^^rj = — 27 ^ ^ 0. 

Theorem 6. If H > then function y(x) = —x is a (local) minimum for 

the unrestricted problem l[T]l. For H < ^r the problem has no solution. 

Remark 7. The unrestricted problem (JQ does not admit global minimum. Take 
indeed, for large values of the parameter a, the control function 



u{x) 



'a if <*<§+§ 



This gives R[u(-)] = 1+ r aS which vanishes as a — > +oo, showing that no global 
solution can exist. 

By the symmetry with respect to the yy axis, a local solution to the unre- 
stricted two-dimensional Newton's problem of minimal resistance with H > r 
is a triangle, with value for resistance R equal to r -i+ H i ■ 

5 Restricted problem 

We now study problem (JTJ with £1 = H.J. In this case the optimal control can 
take values on the boundary of the admissible set of control values fl (u = 0). 
If the optimal control u(') is always taking values in the interior of f2, tt(a;) > 
Vie [0, r], then the optimal solution must satisfy © and it corresponds to the 
one found in 21 

u(x) = — , VarefO.r], (10) 
r 

with resistance 

#=^5 FT?- (11) 

r 2 + H 2 y J 

We show next that this is solution of the restricted problem only for H > r: for 
H < r the minimum value for the resistance is R = r — *t . 

It is clear, from the boundary conditions y(0) = 0, y(r) = H, r > 0, H > 0, 
that u(x) = 0, V i £ [0,r], is not a possibility: there must exist at least one 
non-empty subinterval of [0, r] for which u(x) > (otherwise y(x) would be 
constant, and it would be not possible to satisfy simultaneously y(0) = and 
y(r) = H). The simplest situations are given by 



(12) 




7 



or 

«.)-{.* * "f f' (13) 
10 if £ < ^ < ^ • 

We get (E3J| from JEJ) taking £ = 0; (JUIJ) from JEJl with £ = r. For JEJl the 

resistance is given by = S, + rp^^i^jp , that has a minimum value for 

£ = r - i? > 0: i?(r - ff) = r - f , 

, v JO if 0<.T<r-i/, 
«W = < n ., . (14) 

II if r — H < x < r . 

For r = H (|H|> coincides with (jTU|) ; for r > H 

H\ ( r 3 \ _ H(r — H) 2 
Yj~ \ r 2 + H 2 ) ~ ~2(r 2 +ff 2 ) < °' 

and (|14|l is better than (|10|l . Similarly, for 1(13(1 the resistance is given by 

R (0 = £^ 1 p+r-{, (15) 
that has minimum value for £ = H > 0: 

n(,) = | 1 if (16) 
[0 if ff<a;<r, 

R{H) = r—M-, which coincides with the value for the resistance obtained with 
l)14[l. If one compares directly (|1 l|l with 115(1 one get the conclusion that 1(10(1 
is better than 1(13(1 precisely when r < H: 

r 2 + i? 2 + ^ 2 / [(r - £) 2 + # 2 ] (r 2 + tf 2 ) ' V ' 

and since — iJ 2 < r 2 — r£ — i? 2 < r 2 — i/ 2 , ((T7|l is negative if r < H, that is, for 
r < H JTDJl is better than (JT3J. For r — H (JTHJ) coincide with JTUJ), for r > H 
l(TB|) is better than iJTUjl and as good as l|TI)l . 

We now show that for r > H it is possible to obtain the resistance value 
r — y from infinitely many other ways, but no better (no less value) than this 
quantity. Generic situation is given by 

,s [O if 6i < £ < , i = 0,...,n, 

where n G N, = £ < £i < ■ ■ • < £2^+1 = r, = fi < fi± < ■ ■ ■ < \i n = H. We 
remark that for the simplest case n = 1 ((18(1 simplifies to 

'0 if 0<x<£i, 

"i<-'') = 4 5=5T if & ' •'• ' & • 
if £ 2 < x < r , 
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which covers all the previously considered situations: for £1 = 0, £2 = f we 
obtain (|TU)) : for £2 = r dt i and for £1 = one obtains (JT5J. All Pontryagin 
control extremals of the restricted problem are of the form l|18|) . and by The- 
orem 0] also the minimizing controls. The resistance force R n associated with 
<(T%|) is given by 

Rn (&) ■ ■ • j £2n+l)M0> • • • ) M") 

n n— 1 , t t ,3 

=E(^-f»)+E7> — t +2 ^t — (i9) 

It is a simple exercise of calculus to see that function (11911 has three critical 
points: two of them not admissible, the third one a minimizer. The first critical 
point is defined by /i; = 0, i = 0, . . . , n, which is not admissible given the fact 
that n n = H > 0. The second critical point is given by /ij — fii-x = £21-1 — £,21, 
i = 1, . . . ,n, which is not admissible since fii — fii—i > 0, £21-1 — £2; < 0, and 
jtx, = i = 1, ... ,n, is not a possibility given /x n = H > fio — 0. The third 

critical point is 

Hi - Hi-i = &i - , i = l,...,n, (20) 

which is a minimizer for H < r. Thus, all the minimizing controls for the 
restricted two-dimensional problem with H < r are of the following form: 

, s JO if &t < a; < 5 i = 0,...,n, 

[ 1 if £ 2s :+i < a; < £ 2 i+2 , » = 0, . . . , n - 1 , 

n = 1,2, = £ < 6. < ••• < 6n+i = r< For u n (x) given by (|2TJ> the 
resistance (ftT))l reduces to i?„ = r — Ip, V n G N. 

Theorem 8. The restricted two-dimensional Newton's problem of minimal re- 
sistance admit always a solution: 

• the unique solution associated to control l|10|) . when H > r; 

• infinitely many solutions associated to the controls l|21|l . when H < r. 

3 

In the case H > r the minimum value for the resistance is >± r , H '2 , otherwise 
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6 Conclusion 

Newton's classical problem of minimal resistance offer two interesting situations 
to be studied: the problem in dimension two; and the problem in dimension d, d 
a real number greater or equal than three. While second situation is well studied 
in the literature, and well understood, the first one has been ignored. In the 
classical three-dimensional Newton's problem of minimal resistance, only the 
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problem with restriction u(x) = y(x) > makes sense (without the restriction 
the problem has no local minimum). In the two-dimensional case, we have 
proved that the unrestricted case is also a well defined problem when H > ^-r, 

the minimum value for the resistance being r -i r +H i ■ The local minimizcr is a 
triangle. The two-dimensional problem with restriction u(x) = y(x) > has 
always a solution: a unique solution (a triangle) when H > r, with value for 
resistance equal to the unrestricted case; infinitely many alternative solutions 
for r > H, the minimal aerodynamical resistance being r — ^. 
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